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We study the stochastic dynamics of c and b quarks, produced in hard initial processes, 
in the hot medium created after the collision of two relativistic heavy ions. This is done 
through the numerical solution of the relativistic Langevin equation. The latter requires 
the knowledge of the friction and diffusion coefficients, whose microscopic evaluation is 
performed treating separately the contribution of soft and hard collisions. The evolution 
of the background medium is described by ideal/viscous hydrodynamics. Below the 
critical temperature the heavy quarks are converted into hadrons, whose semileptonic 
decays provide single-electron spectra to be compared with the current experimental 
data measured at RHIC. We focus on the nuclear modification factor Raa an d on the 
elliptic-flow coefficient V2, getting, for sufficiently large px, a reasonable agreement. 



1. Introduction 

Heavy quarks, produced in initial hard processes, allow to perform a "tomography" 
of the medium (hopefully a Quark Gluon Plasma) created in heavy-ion collisions. 
The modification of their spectra provides information on the properties (encoded 
into few transport coefficients) of the matter crossed before hadronizing and giving 
rise to experimental signals (so far, the electrons from their semi-leptonic decays). 
We employ an approach based on the relativistic Langevin equation, assuming 
that medium-modifications of the initial heavy-quark spectrum (in particular its 
quenching at high px) arise from the cumulative effect of many random collisions. 
The final goal of our analysis is to achieve a reasonable description of the Raa and V2 
of non-photonic electrons from heavy- flavor decays measured at RHIC. Preliminary 
results can be found in Ref. For similar studies see also Refs. I 2PK I 5 I 



2. The Langevin equation in a dynamical medium 

We solve the Langevin equation for heavy quarks propagating in the expanding 
fireball created in heavy-ion collisions. We assume that local thermal equilibrium is 
reached, so that the background medium is entirely described by the four-v elocity 
and temperature fields ^(x) and T(x), provided by two different hydro codes ISLlIsI 
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This requires a generalization of the algorithm used in Ref. ^ for the static case. 
We focus on a given quark at (x n ,p n ) after n steps of evolution. Next we move to 
the local fluid rest-frame and update its position and momentum by the quantities 
Ax n = (p n /E p )Ai and 

A#, - -r, D (p n )f n M + e(i)At = -n D {p n )f n M + 9 v(p n )(i(t)VAl, (1) 

where we take - in the fluid rest-frame - At = 0.02 fm/c. In the above we express 
the noise term through the tensor (we omit the "bar" ) 

g» (p) = y/^pjp 4 ? + y/^pjtf* - ffP ) , (2) 

depending on the transverse/longitudinal momentum diffusion coefficients kt/l{p)) 
and the uncorrelated random variables with (C(t)C 3 it')) =S l ^5tt'. Hence, one 
simply needs to extract three random numbers £ J from a gaussian distribution with 
cr = l. Then, going back to the Lab frame, one obtains the updated (x n+ i,p n+1 ). 



3. Evaluation of the transport coefficients 

The coefficients 

yield the average transverse/longitudinal squared- momentum acquired pe r u nit 
time by the heavy quark through the collisions in the medium. Following Ref. ^3 we 
introduce an intermediate cutoff \t\* ^m 2 D (t=(P' ' — P) 2 ) separating hard and soft 
scatterings. The contribution of hard collisions (\t \ > \t\*) is evaluated through a ki- 
netic pQCD calculation of the processes Q(P)q i /i—^Q(P')q i /i and Q(P)g—>Q(P')g. 
On the other hand in soft collisions (\t\ < \t\*) the exchanged gluon has "time" to 
feel the presence of the other particles. A resummation of medium effects is thus 
required and this is provided by the Hard Thermal Loop approximation. The final 
result is given by the sum of the two contributions n T / L {p) — Ky^(p) + n^^ip) 

and its explicit expression can be found in Ref. In Fig. [T] we display its behavior 
in the case of a c quark in a plasma at T = 400 MeV. The sensitivity to the value 
of the intermediate cutoff \t\* is quite small, hence supporting the validity of the 
approach. Finally, as discussed in Ref. ®, the friction coefficient r)u{p) is fixed so 
that the heavy quarks evolve towards thermal equilibrium. 



4. Numerical results 

For each explored case we generated an initial sample of 45 • 10 6 cc and bb pairs, 
using the POWHEG code with CTEQ6M PDFs. In the AA case we intro- 
duced nuclear effects in the PDFs according to the EPS09 scheme the quarks 
were then distributed in the transverse plane according to the nuclear overlap 
function dN/dx±^TAB(x,y) = Ta{x + b/2, ij)Tb{x — 6/2, y) and given a further 
kx broadening on top of the "intrinsic" one. At the proper-time r = t 2 — z 2 = tq 
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Fig. 1. The momentum-diffusion coefficients k t ^ l of a c quark after summing the soft and hard 
contributions. The dependence on the cutoff |t|* is mild. The coupling g was evaluated at /i = 2ttT. 





CTcc (fJ-b) 




Hydro code 


T (fm/c) 


so (fm- 3 ) 


T (MeV) 


PP 


254.14 


1.769 


ideal 


0.6 


110 


357 


AA 


236.11 


2.033 


viscous 


1.0 


83.8 


333 



Table 1. Initialization for RHIC: the cc and bb production cross section at v^jvjv = 200 GeV given 
by POWHEG and the explored hydro scenarios. 



we started following the Langevin dynamics of the quarks unt il hadronization. The 
latter was modeled using Peterson fragmentation f unction s with branching frac- 
tions into the different hadrons taken from Refs. | 15 | 16 [ Finally each hadron was 
forced to decay into electrons with PYTHIA ^3 using updated decay tables ^Sl. 
The e-spectra from c and b were then combined with a weight given by the re- 
spective total production cross-section quoted in Table Q] where the hydro initial 
conditions are also reported. We first consider the nuclear modification factor 
Raa{pt) = {dN/d P T) AA /(N co ii)(dN/dp T ) pp of the heavy-flavor spectra. In Fig. H 
we display the effects of fragmentation and semi-leptonic decays on the charm Raa, 
both leading to an additional quenching of the spectrum. We also show the sepa- 
rate contribution of the electrons from charm and bottom. In Fig. [3] the Raa of 
non-photonic electrons measured by Phenix is compared with our findings for dif- 
ferent values of the coupling. The high-momentum region (pt~4 GeV/c) is better 
reproduced evaluating a s at the scale 1.57rT. The dependence on the hydro scenario 
(ideal/viscous) turns out to be small. In Figs. 0] and [5] the same kind of analysis 
is performed for the elliptic-flow coefficient i>2 (pt) = (cos(20))p r , with 4> measured 
with respect to the reaction plane. Our results for the electron w 2 underestimate 
the experimental data. The same occurs for Raa at moderate pr- Part of this dis- 
agreement could arise from ignoring coalescence as a mechanism of hadronization, 
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Fig. 2. Left panel: the effect of fragmentation and semi-leptonic decays on the charm -RaAj both 
leading to a quenching of the spectrum. Right panel: the separate contribution of charm and 
bottom to the single-electron Raa- Electrons from fo-decays start contributing for py?52 GeV/c. 




Fig. 3. The electron Raa for minimum bias Au-Au collisions at = 200 GeV. The sensitivity 

to the scale fi = nT — 2irT is displayed. The high- momentum region (py ^ 4 GeV/c) is better 
reproduced with the intermediate coupling. Ideal (left panel) and viscous (right panel) hydro 
scenarios provide similar results. 



as can be inferred from the analysis of its effect performed in Ref. ^ A possible 
further discrepancy would require a better evaluation of the transport coefficients. 
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Fig. 4. Left panel: the effect of fragmentation and semi-leptonic decays on the elliptic-flow of 
charm, resulting in a shift to lower values of px- Right panel: the separate contribution of charm 
and bottom to the single-electron i>2- The bottom does not display any flow. 
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Fig. 5. The elliptic-flow coefficient V2 of non-photonic electrons for minimum bias Au-Au collisions 
at y / i NN = 200 GeV. The sensitivity to the scale /i = irT — 2itT is displayed. The data would favor 
the strongest value of the coupling, which however would provide a too small value of Raa- Ideal 
(left panel) and viscous (right panel) hydro scenarios provide again similar results. 
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